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Abstract 

A  new  design  has  recently  been  proposed  in  the  field  of  solid  oxide  fuel  cells,  consisting  of  a  traditional  electrochemical  cell  integrated 
with  a  pre-heater.  In  this  paper  a  simulation  model  for  the  rectangular  planar  solid  oxide  fuel  cell  with  integrated  air  pre-heater  is  presented. 
A  two-dimensional  stack  simulation  is  presented  as  well,  one  axis  coincides  with  the  fuel  fiow  direction,  the  other  with  the  stack  height.  Local 
quantities  such  as  current  density,  gas  and  solid  temperatures  are  reported  and  cell  characteristics  predicted.  In  a  parameter  study,  effects  of 
oxygen  utilisation  and  heat-transfer  conditions  in  the  pre-heater  on  the  local  temperature  distribution  of  the  solid  structure  are  considered.  As 
a  result,  the  benefit  of  the  new  cell  design  becomes  evident  when  low  air  flow  rates  are  applied.  A  further  advantage  associated  with  the 
reduced  flow  rate  is  the  low  air  temperature  at  the  inlet.  ©  1997  Elsevier  Science  S.A. 
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1.  Introduction 

Ceramic  materials  of  a  solid  oxide  fuel  cell  (SOFC)  [  1  ], 
i.e.  anode,  cathode  and  electrolyte,  cannot  withstand  steep 
temperature  gradients.  Thus,  in  order  to  obtain  an  almost 
uniform  temperature  distribution  on  the  cer'  plane,  a  higher 
air-flow  rate  than  that  required  for  fuel  oxidation  must  be  fed 
to  the  cell  at  high  temperature.  As  a  consequence  the  air  pre¬ 
heater  and  the  air  blower  become  the  most  expensive  com¬ 
ponents  of  the  whole  fuel  cell  plant:  for  example,  they  account 
for  about  65%  to  the  investment  cost  [2]  when  external 
methane-steam  reforming  is  applied. 

In  order  to  overcome  this  problem  a  rectangular  SOFC 
with  integrated  air  pre-heater  has  recently  been  proposed  by 
KFA  Jiilich  [3].  A  circular  shape  for  the  cell  and  the  stack 
with  integrated  air  pre-heater  has  been  proposed  by  others 
[4J,  and  will  not  be  considered  here. 

The  basic  idea  of  this  design  is  to  shift  the  temperature 
gradient  Tout  -  Tin  from  the  air  channels  of  a  conventional  cell 
to  the  air  flowing  into  the  integrated  heat  exchanger.  Thus, 
the  air-inlet  temperature  could  be  reduced  by  the  same 
amount,  or  even  more  if  the  air-mass  flow  is  decreased. 

The  heat-exchanger  integrated  stack  is  formed  by  a  series 
of  repeat  modules,  each  of  them  constituted  of  a  traditional 
planar  electrochemical  cell  coupled  with  a  heat-exchanging 
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element.  A  detailed  design  of  the  faditional  planar  cell,  show¬ 
ing  more  details  such  as  the  gas  channels  that  distribute  fuel 
and  oxidant  on  the  cell,  can  be  found  elsewhere  [  1  ] .  A  sketch 
of  a  rectangular  heat-exchanging  element  is  shown  in 
Fig.  1(a):  it  is  a  void  chamber  where  the  air  flows  before 
entering  the  cathodic  compartment  of  the  electrochemical 
cell.  No  channels  are  present.  Finally,  a  scheme  of  the  rec¬ 
tangular  stack  is  presented  in  Fig.  1(b). 

A  further  proposal  of  innovation  that  will  be  discussed  in 
the  present  paper  is  the  use  of  metallic  materials  instead  of 
ceramic  components  for  the  interconnector  and  for  the  inte¬ 
grated  heat  exchanger. 

The  fuel  (pure  hydrogen)  and  the  cathodic  air  stream  are 
in  co-flow,  whereas  the  air  is  pre-heated  in  counter-flow  to 
the  cathodic  air.  Before  entering  the  fuel  cell  fresh  air  passes 
into  the  afterburner  compartment.  Here,  the  fuel  not  con¬ 
verted  in  the  cell  is  burnt  to  contribute  to  the  air  pre-heating. 
In  this  way  the  fuel  utilisation  is  closely  related  to  the  air 
excess  and  therefore  subject  to  optimisation  procedures. 

Previous  types  of  SOFC  simulation  models  [  5— 17]  have 
to  be  extended  to  account  for  the  additional  internal  heat- 
transfer  process  in  the  integrated  heat  exchanger.  In  the  pres¬ 
ent  work  a  model  accounting  for  the  presence  of  the  integrated 
heat  exchanger  is  discussed  and  on  this  basis  a  new  simulation 
code  is  presented  in  application  to  a  rectangular  single  cell 
and  to  a  stack  as  well.  Computational  results  are  discussed 
and  the  benefit  of  this  new  design  is  pointed  out. 
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I;ig.  I .  { a >  View  of  the  integrated  heat-exchanger  plate.  ( b )  Drawing  of  the  stack  1 3 1 . 


2.  Simulation  model 

2. 1 .  Hypotheses  of  the  model  and  input  data 

The  present  stack  simulation  program  is  two-dimensional 
in  space  and  is  time-dependent.  The  length  co-ordinate,  .v,  is 
parallel  to  the  fuel  stream  direction,  the  origin  coincident  with 
the  point  where  the  fuel  enters  the  cell.  The  ^-coordinate  is 
perpendicular  to  the  cell  plane.  A  one-dimensional  model  of 
stacked  cells  is  allowed  because  co-  or  counter-flow  occurs. 
This  implies,  however,  that  the  outer  stack  surfaces  in  the  a- 
-  plane  must  be  considered  adiabatic.  The  remaining  surfaces 
are  allowed  to  radiate  or  convey  heat  to  the  surrounding  shell 
which  is  assumed  to  have  a  uniform  and  constant  temperature. 

As  the  current  collectors  have  a  very  high  electric  conduc¬ 
tivity.  they  can  be  considered  as  equipotential  plates. 

The  following  quantities  are  required  as  input  data: 

1 .  Design  parameters :  ( i )  size  and  configuration  of  the  elec¬ 
trolyte,  electrodes,  interconnector,  heat-exchanger  plate: 
(ii)  size  and  number  of  gas  channels;  (iii)  number  of 
cells. 

2.  Operating  conditions:  (i)  pressure,  inlet  temperature, 
composition  and  utilisation  of  the  anodic  and  cathodic 
gases;  ( ii )  shell  temperature;  ( iii )  average  current  density 
or  cell  voltage, 

3.  Properties  and  empirical  data:  ( i )  definition  of  the  mate¬ 
rials  used;  (ii)  empirical  parameters  for  the  rate-limiting 
kinetics. 

The  main  output  quantities  are:  ( 1 )  electrical  power  and 
efficiency;  (2)  local  distributions  of  current  density,  gases 
and  solid  temperatures,  gas  compositions  and  flowrates. 


2.2.  Balance  equations 

2.2.  /.  Simplifications 

As  already  discussed  in  Refs.  [  10,16]  the  heat  and  mass 
transfer  in  the  gas  channels  can  be  treated  as  being  one¬ 
dimensional  without  undergoing  considerable  errors.  Fur¬ 
thermore,  the  streamwise  diffusion  terms  can  be  neglected 
since  the  convective  mass  and  heat  fluxes  are  the  prevailing 
contributions.  Thus,  the  mass  and  energy  balances  of  the  gas 
phase  can  be  written  in  plug-flow  form. 

The  streamwise  molar  flux  gradient  is  merely  depending 
on  the  conversion  rates.  Similarly  the  streamwise  temperature 
distribution  is  determined  by  the  convective  heat  transfer 
trom  the  channel  walls,  while  the  solid  temperature  depends 
on  the  reaction  enthalpies  and  on  the  convective  heat 
exchange  with  the  gases. 

2.2.2.  Mass  balance 

Under  the  assumptions  made,  the  mass  balance  can  be 
written  as 

(t/i,  _ 

a7_#'  (»> 

where  i  denotes  the  particular  gas  component  converted. 

The  reaction  rate  can  be  calculated  by  means  of  Faraday’s 
law 

I 


where  ne  is  the  number  of  electrons  transferred  in  the  elec¬ 
trochemical  reaction  and  v  is  the  stoichiometric  factor  of  the 
reacting  component;  v  is  negative  for  reactants  and  positive 
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for  products.  In  the  case  of  a  cell  fed  with  hydrogen  ne  is 
equal  to  2,  corresponding  to  the  following  reactions: 

1  /20z  +  2e  ~  -» O2 "  cathodic  reaction 

H2  +  O2  “  -»  H20  +  2e  ~  anodic  reaction 

Hi  +  1  /20i  -*  HiO  over-all  reaction 


2.2.3.  Energy  balance  of  the  gaseous  phase 
As  mentioned  above,  the  streamwise  temperature  gradient 
in  the  gas  flows  is  assumed  to  depend  only  on  the  convection 
phenomenon  due  to  both  the  heat  and  the  mass  transfer  from 
the  channel  walls  to  the  bulk  of  gases 


E'Vn/jf  +  X>r.§<  T,-T,)  +  heB„k  r,-T.)-0  ( 3 ) 


ojc 


where  B  is  the  ratio  between  the  gas-solid  heat-exchange  area 
and  the  cell  area,  h  the  heat-transfer  coefficient  calculated 
from  the  Nusselt  number,  Nu 


h  =  Nurr  (4) 

« h 

The  value  of  the  Nusselt  number  generally  depends  on  the 
channel  geometry,  the  Reynolds  number.  Re ,  and  the  Prandtl 
number,  Pr.  For  the  laminar  established  flow,  however,  which 
holds  for  the  fuel  cell  channels  the  Nusselt  number  is  a 
constant. 

Eq.  (3)  is  also  applied  to  the  flow  through  the  heat- 
exchanger  section.  In  this  case,  of  course,  the  conversion  term 
vanishes. 


2.2.4.  Energy  balance  of  the  solid  structure 
The  energy  balance  of  the  solid  material  describes  the 
unsteady  heat  conduction  in  the  quasi-homogeneous  solid 
structure  of  the  stack.  The  convective  heat,  transferred  from 
the  gas  flow  to  the  solid,  and  the  reaction  enthalpies  occur  as 
source  terms 


a2rs  „  a ts 

+*=9r 


(5) 


where  Kx  and  K.  are  the  effective  thermal  conductivities  of 
the  solid  structure,  taking  the  particular  cell  configurations  in 
both  directions  into  account.  In  the  ^-direction  parallel  heat 
conduction  through  the  h  parallel  layers  is  assumed 


h 


(6) 


while  for  the  ^-direction  the  heat  conduction  is  in  series 


Kz= 


is* 


h 


(7) 


In  both  effective  conductivities  the  contributions  of  the  gas 
layers  are  also  included. 

2.3.  Local  reaction  rate 


The  internal  electric  resistance  of  the  cell  has  been  deter¬ 
mined  according  to  the  procedure  proposed  by  Refs. 
[5,10,17,18]. 

In  particular,  the  open-circuit  voltage  V™  is  given  by 
Nemst’s  equation  applied  to  the  over-all  reaction 


y°c 


-AG 

neF 


-  AG°  RTS 
neF  neF 


In(npr) 

i 


(8) 


Ionic  resistance  of  the  electrolyte,  ohmic  losses  of  the  elec¬ 
trodes  and  of  the  bipolar  plate  and  finally  contact  resistances 
have  been  considered.  All  quantities  except  contact  resis¬ 
tances  are  expressed  by  temperature  dependent  equation: 


(9) 


Concentration  polarization  is  usually  negligible  for  SOFC 
devices,  while  activation  polarization  at  the  anode  and  cath¬ 
ode  electrodes  has  been  considered,  according  to  the 
equations: 


(10) 

(11) 


So,  the  potential  between  the  current  collectors  is  given  by 


y  cc  _  y  oe _ y  ohm _  y  P°* 


(12) 


where  Vcc  is  the  voltage  at  the  current  collectors,  V/ohl"  and 
V1™1  are  the  voltage  drops  due  to  ohmic  and  polarisation 
effects,  respectively. 


2.4.  Boundary  conditions 

As  shown  in  Fig.  1  the  exhaust  anode  gas  is  burnt  together 
with  the  cathodic  off-gas  in  a  separate  compartment  contain¬ 
ing  the  after-burner  heat  exchanger  which  pre-heats  the  ambi¬ 
ent  air  up  to  the  inlet  temperature  of  the  integrated  pre-heater. 
Heat  is  allowed  to  be  exchanged  by  convection  and  radiation 
from  the  stack  to  this  region.  As  a  simple  assumption  the 
temperature  of  this  compartment  is  considered  to  be  uniform. 
The  heat  transferred  to  the  stack  by  convection  is  given  by 

^  =  a(7’shcii-7’s) 

The  radiative  contribution  is  approximated  by 

v=4efLi^L^  (7^1, -7-,)  (14) 

Because  the  heat-transfer  problem  is  treated  one-dimension- 
ally  in  each  z-plane,  the  side  walls  must  be  considered 
adiabatic. 
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2.5.  Numerical  integration 


3.  Computational  results 


The  balance  Eqs.  ( I ),  (3)  and  (5)  have  to  be  solved  for 
two  spatial  coordinates  and  one  temporal  one.  In  fact,  how¬ 
ever,  it  is  sufficient  to  treat  only  Eq.  (5)  as  time-dependent, 
since  mass  transfer  and  all  chemical  and  fluid  dynamics  proc¬ 
esses  are  fast  compared  with  the  heat  dispersion  in  the  solid 
structure.  Thus,  a  steady -state  solution  of  the  mass  and  energy 
balances  of  the  gaseous  phase  is  found  at  each  time  step. 
Successively  a  new  solid-temperature  distribution  is  calcu¬ 
lated  using  the  source  terms  of  the  preceding  time  step. 

The  integration  is  performed  by  applying  finite  differenc¬ 
ing.  Fast  convergence  for  the  solution  of  Eq.  ( 5 )  is  achieved 
by  means  of  the  successive  overrelaxation  method  (SOR). 
The  corresponding  SOR  factor  was  automatically  adjusted 
during  the  running  of  the  program,  proportionally  to  the  dis¬ 
tance  from  the  final  solution. 

In  particular,  the  following  procedure  was  applied  to  solve 
the  equations  for  the  single  cell: 

1.  specify  the  operating  conditions  in  the  input  file  and  set 
an  initial  solid-temperature  distribution; 

2.  calculate  the  current  density  and  solve  the  mass  balance; 

3.  solve  the  energy  balance  of  the  gaseous  phase  to  determine 
the  temperature  distribution  along  the  gas  channels,  and 

4.  solve  the  energy  balance  of  the  solid  to  find  a  new  tem¬ 
perature  field. 

Steps  I  through  4  are  repeated  until  the  steady-state  solution 
is  reached  within  a  given  threshold  of  accuracy. 

The  calculation  of  the  whole  stack  starts  from  the  solution 
of  the  single  cell.  The  computation  proceeds  as  follows: 

1 .  extend  the  results  of  the  single  cell  to  the  whole  stack; 

2.  solve  mass  and  energy  balances  of  the  gas  flow  tor  each 
cell,  and 

3.  solve  the  energy  balance  of  the  solid,  accounting  for  the 
outer  boundary  conditions. 

Repeat  steps  2  and  3  until  the  steady-state  solution  is  reached 
within  the  accuracy  required. 

The  average  temperature  level  of  the  solid  material  is  an 
important  design  quantity.  Since  it  strongly  depends  on  the 
inlet  temperature  and  mass  flow  of  the  air,  an  iteration  step 
seems  to  be  useful  to  find  out  the  solution  satisfying  the 
average  required  value  of  the  solid  temperature.  For  this  pur¬ 
pose  a  PID  controller  ( proportional,  integral,  differential )  for 
the  air-inlet  temperature  is  simulated.  Starting  from  an  esti¬ 
mated  temperature  field  of  the  solid,  the  temperature  of  the 
inlet  air  is  varied  proportionally,  integrally  and  differentially 
to  the  difference  between  the  actual  average  temperature  of 
the  solid  and  the  target  value. 

The  computer  code  has  been  applied  to  a  single  adiabatic 
cell  and  a  stack  consisting  of  five  cells.  For  each  computation 
presented  in  this  paper  an  average  solid  temperature  of 
Tav  =  1 193  K  was  required  which  results  in  different  air-inlet 
temperatures  at  the  integrated  heat  exchanger. 


3. 1.  The  single  adiabatic  cell 

Figs.  2  and  3  show  the  temperature  distributions  of  the 
fuel,  of  the  air  and  of  the  solid  structure  along  a  cell  operating 
at  an  average  current  density  of /= 3000  Am-2.  The  fuel  is 
hydrogen,  80%  of  which  is  utilised.  The  stoichiometry  factor 
of  the  air  flow  is  A  =  2.7  resulting  in  an  02  utilisation  of 
Uo,  =  0.3.  Further  design  and  operating  data  are  listed  in 
Table  I. 

Fig.  2  represents  the  temperature  distribution  of  the  fuel, 
of  the  cathodic  air  and  of  the  air  in  the  integrated  heat 
exchanger  along  the  cell.  To  come  out  with  an  average  solid 
temperature  of  7s  av  =1193  K  the  air  has  to  enter  the  stack 
with  a  temperature  of  jTair  j„  =  7 10  K.  The  air  stream  running 
in  counter-flow  to  the  cathodic  air  is  continuously  heated  up 
by  the  waste  heat  electrochemical  ly  produced.  The  tempera¬ 
ture  of  the  cathodic  and  the  anodic  gases  are  close  to  that  of 
the  solid  material. 

Fig.  3  exhibits  the  solid-structure  temperature  in  a  more 
expanded  scale.  The  maximum  temperature  is  about 


Fig,  2.  Temperatures  ot  the  fuel,  the  cathodic  air  and  the  air  flowing  into  the 
integrated  heat-exchanger  plate  along  the  cell  length. 


Fig.  3.  Distribution  of  temperature  on  the  solid  as  a  function  of  the  cell 
length. 
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Table  I 

Design  parameters  and  operating  conditions  for  the  simulated  rectangular 
single  cell 


Simulated  rectangular  cell 

Design  parameters 

Overall  cell  dimensions 

0.12  mXO.l  m 

Active  cell  area 

0.01  nr 

Anode,  cathode,  electrolyte  thickness 

1  X  I0“4  m 

Heat-exchanger  plate  thickness 

5X10  ’m 

Interconnector  thickness 

2.4X10  'm 

Anodic /cathodic  channel  shape 

square 

Channel  dimension 

1  X  10“ -'m 

Rib  dimension 

1  X  10  '  m 

Operating  conditions 

Current  density 

3000  Am- 

Cell  voltage 

0.71  V 

Pressure 

1.013  X  10  s  Pa 

Anodic  inlet  temperature 

1 173  K 

Cathodic  inlet  temperature 

711  K 

Average  solid  temperature 

1 193  K 

Anodic  inlet  composition 

97%H,  +  3f/rH,0 

Cathodic  inlet  composition 

air 

Hydrogen-utilisation  factor 

0.8 

Oxygen-utilisation  factor 

0.3 

Heat  exchange 

adiabatic  conditions 

Ts  —  1220  K  and  the  difference  observed  across  the  cell  is 
about  A7>70K  which  is  a  small  value  with  a  view  to  the 
low  air  stoichiometry  factor  of  A  =  2.7. 

Due  to  the  low  mass  flux  the  high  thermal  conductivity  of 
the  solid  material  contributes  considerably  to  removal  of  heat 
from  the  entrance  region  of  the  cell  to  the  exit  where  the  heat 
is  released  to  the  refrigerant  air.  Then  the  heat  removal  from 
areas  of  high  waste  heat  production  is  not  predominantly  done 
by  convection,  while  the  heat  conduction  through  the  solid 
structure  becomes  the  prevailing  mechanism. 

The  local  current-density  distribution  is  displayed  in  Fig.  4. 
As  a  result  of  the  high  fuel  utilisation  the  data  vary  in  a  wide 
range  from  about  7000  A  m  ~ 2  at  the  entrance  to  1 200  Am-2 
at  the  outlet.  With  decreasing  utilisation,  of  course,  the  uni¬ 
formity  would  improve. 


cell  length  [m] 

Fig.  4.  Electrical  current  density  distribution  along  the  cell  length. 


Fig.  5.  Characteristic  curve  for  the  single  cell. 


In  Fig.  5  the  characteristic  curve  of  the  adiabatic  single  cell 
is  reported.  The  average  cell  temperature  is  kept  constant,  but 
different  solid  temperature  distributions  arise  at  different 
operating  conditions.  Small  variations  of  the  cell  resistance 
take  place  and  cause  small  changes  of  the  line  slope,  which, 
anyway,  is  very  similar  to  a  typical  SOFC  straight  line,  the 
voltage  of  the  cell  ranging  from  0.71  V  at  3000  A  m'2  to 
0.81  V  at  1000  Am-2. 

In  a  conventional  cell  built  with  ceramic  materials  and  not 
equipped  with  an  integrated  heat  exchanger  maximum  solid- 
and  gas-outlet  temperatures  of  Ts  =  1 370  K  would  occur,  even 
if  the  inlet  air  temperature  is  reduced  to  985  K.  The  corre¬ 
sponding  minimum  solid  temperature  is  Ts=1060  K,  as 
shown  in  Fig.  6.  The  improvements  due  to  the  use  of  metallic 
interconnectors  are  pointed  out. 

Fig.  7  shows  how  the  heat-transfer  conditions  in  the  inte¬ 
grated  heat  exchanger  affect  the  solid  temperature  distribu¬ 
tion.  If  the  label  ‘Tin  Fig.  7  stands  for  the  heat-transfer 
coefficient  calculated  on  the  basis  of  the  data  reported  in 
Table  1,  the  remainder  labels  indicate  by  which  fraction  this 
value  is  diminished  in  the  simulation.  Thus  the  response  of 
the  cell  solid  temperature  to  the  efficiency  of  the  integrated 


Fig.  6.  Comparison  between  the  solid-temperature  distributions  obtained  for 
different  cell  configurations:  with  or  w-thout  integrated  heat  exchanger,  the 
latter  simulated  with  both  a  ceramic  and  a  metallic  interconnector.  Geomet¬ 
rical  and  operating  parameters  are  the  ones  reported  in  Table  1  in  all  cases. 
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Fig.  7.  Effect  of  variations  of  the  heat-exchange  coefficient  for  the  air  flowing 
into  the  integrated  exchanger  on  the  solid-temperature  distribution. 


Tuble  2 

Design  parameters  and  operating  conditions  tor  the  simulated  rectangular 
stuck 


Simulated  rectangular  stack 

Design  parameters 

Nos.  of  cells 

5 

Operating  conditions 

Current  density 

3000 Am  ’ 

Stack  voltage 

3.55  V 

Pressure 

1.013  x  10*  Pa 

Anodic  inlet  temperature 

1 173  K 

Cathodic  inlet  temperature 

191  K 

Average  solid  temperature 

II93K 

Shell  tcmperatuic 

1 183  K 

Anodic  inlet  composition 

‘>7m+3mo 

Cathodic  inlet  composition 

air 

Hydrogen-utilisation  factor 

0,8 

Oxygen-utilisation  factor 

0.3 

Heat  exchange 

radiation + convection 

heat  exchanger  is  illustrated.  For  increasing  label  numbers 
the  operating  conditions  come  continuously  closer  to  those 
of  a  cell  without  integrated  heat  exchanger. 

If  the  heat-transfer  coefficient  is  reduced  by  a  factor  of  8 
the  temperature  distribution  of  the  solid  structure  is  almost 
flat.  By  this  measure  the  pre-heating  air  is  more  slowly  heated 
up  and  thus  able  to  remove  the  large  amount  of  waste  heat  in 
the  entrance  region  of  the  electrochemically  active  part  of  the 
cell,  where  heat  is  mostly  dissipated. 

A  reduction  of  the  heat-transfer  coefficient,  for  instance, 
can  easily  be  verified  by  increasing  the  channel  height  of  the 
integrated  heat  exchanger. 

As  a  matter  of  fact  this  figure  does  point  out  the  benefit  of 
the  integrated  heat  exchanger:  it  is  sufficient  to  optimise  the 
design  of  the  integrated  heat-exchanger  plate  to  obtain  a  flat 


temperature  profile.  This  optimisation  depends  on  the  oper¬ 
ating  conditions  (i.e.  the  02-utilisation factor). 

In  the  case  of  a  traditional  cell  ( i.e.  without  integrated  heat 
exchanger)  the  only  way  to  achieve  such  a  result  would 
consist  of  an  increase  in  the  dimensions  of  the  gas  channels 
that  distribute  fuel  and  oxidant  to  the  cell  plane.  In  this  way 
the  path  of  the  electrons  through  the  electrodes  and  the  dif¬ 
fusion  of  gases  would  become  more  and  more  difficult,  and 
the  electrical  resistance  of  the  cell  would  increase.  For  exam¬ 
ple,  if  the  rib  side  and  the  channel  side  lengths  were  2  x  10" 3 
m  instead  of  1  X  10“ 3  m  (as  it  is  assumed  in  the  model. 
Table  1 ) .  diffusion  problems  of  the  hydrogen  in  the  electrode 
would  take  place:  the  PEN  electrical  resistance  would 
become  about  three  times  greater  [19J  (PEN:  positive  + 
electrolyte  +  negative,  i.e.,  anode  +  electrolyte  +  cathode ) . 

As  the  interconnector  and  the  integrated  heat  exchanger 
have  an  electrical  resistance  that  is  negligible  with  respect  to 
the  PEN  one,  any  change  in  their  geometry  would  have  neg¬ 
ligible  consequences  on  the  electrochemical  performances  of 
the  cell.  The  optimisation  of  integrated  heat-exchanger  cells 
would  be  recommended,  also  for  cells  using  only  ceramic 
materials. 

To  show  the  advantages  of  the  new  design,  and  in  particular 
the  benefit  obtained  by  means  of  the  optimisation,  a  parameter 
study  was  made  investigating  the  influence  of  the  02-utilis¬ 
ation  on  the  solid  temperature  and  on  the  air  temperature  at 
the  inlet.  The  results  for  the  heat-exchanger  integrated  cell 
and  the  conventional  cell  are  plotted  in  Figs.  8  and  9.  The 
new  design  has  been  optimised  for  this  comparison  with 
respect  to  the  heat-exchanger  conditions  assuming  the  heat- 
exchange  coefficient  to  be  devided  by  8. 

The  figures  show  that  the  temperature  variations  along  the 
cell  increase,  while  the  inlet  temperature  decreases,  with 
increasing  02-  utilisation.  As  long  as  the  temperature  profile 
of  the  cell  becomes  flatter,  the  air-outlet  temperature,  and  as 


a  consequence  the  inlet  one,  increases  with  respect  to  the  non 
optimised  case  (Fig.  9). 

•100 

- traditional  cel 

. integrated  heat  exchanger  cel  -  non  optirrized 

- integrated  he*  exchanger  cel  -  optirrized 
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Fig.  8.  Effect  of  variations  of  the  iniet  air  flow  rate  on  the  solid-temperature 
gradient  along  the  cell. 
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UoaH 

Fig.  9.  Effect  of  variations  of  the  inlet  air  flow  rate  on  the  inlet  air  temperature 
necessary  to  keep  the  average  temperature  of  the  cell  at  1 193  K. 


3.2.  Five-cell  snick 

The  preceding  parameter  studies  have  been  performed  for 
a  single  cell  which  may  be  regarded  as  a  centrally  positioned 
cell  in  a  large  stack.  In  order  to  show  the  effect  of  the  outer 
boundary  conditions  on  the  stack  operation  a  stack  consisting 
of  five  cells  is  considered.  Details  of  the  operating  conditions 
are  summarised  in  Table  2.  It  is  assumed  that  the  end  plates 
radiate  and  convey  heat  to  the  surroundings  kept  at  a  tem¬ 
perature  of  rslM!, |  —  1 1 83  K. 

Fig.  10  shows  the  temperature  distribution  of  the  solid 
structure  along  the  five  cells  of  the  stack.  The  temperature 
level  of  the  end  cells  is  lower  than  that  of  the  central  cells. 
Since  for  cell  1  —  in  the  drawing  of  Fig.  1  the  left  one  —  the 
outer  channel  is  an  air  channel  and  for  cell  5  a  fuel  channel, 
the  temperature  distributions  of  the  end  plates  are  different. 
The  higher  mass  flow  on  the  air  side  supplies  cell  1  with  heat, 
so  that  the  cell  temperature  does  not  drop  so  much  —  com¬ 
pared  with  cell  5  —  by  the  effect  of  radiative  heat  removal 
from  the  ceil  surface. 


Fig.  10.  Solid-temperature  distribution  along  the  cell  co-ordinate  for  the 
stack. 


Fig.  1 1 .  Electrical  current-density  distribution  along  the  cell  co-ordinate  for 
the  various  cells  of  the  slack. 


Fig.  1 2.  Characteristic  curve  for  the  rectangular  stack. 


Due  to  the  small  temperature  variations  across  the  stack, 
only  minor  variations  of  the  local  current  density  occur 
between  the  particular  cells  as  demonstrated  in  Fig.  1 1 . 

Finally,  Fig.  12  represents  the  stack  characteristic.  For 
j  =  3000  A  m  ~ 2  a  stack  voltage  of  about  U  =  3.6  V  is  predicted 
corresponding  to  an  average  cell  voltage  of  (7=0.71  V.  The 
variation  of  this  quantity  along  the  stack  is  not  significant. 


4.  Conclusions 

The  operating  conditions  of  a  rectangular  planar  SOFC 
with  integrated  air  pre-heater  have  been  analysed.  The  idea 
of  this  new  fuel  cell  type  is  to  reduce  the  cost  of  the  expensive 
external  air  pre-heater  by  applying  lower  air  flow  rales 
together  with  lower  inlet  temperatures. 

The  parameter  study  indicates  that  the  benefit  increases  for 
low  air-excess  rates,  i.e.  for  high  oxygen  utilisation.  For  this 
case  the  air-inlet  temperatures  can  drastically  be  reduced.  The 
air  in  the  pre-heater  streams  in  counter-flow  to  the  cathodic 
air  and  is  heated  up  by  the  electrochemically  produced  waste 
heat.  The  high  thermal  conductivity  of  the  metallic  solid 
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structure  contributes  to  the  compensation  of  the  temperature 
gradients  occurring  along  the  cell.  Thus,  the  heat  produced 
in  the  entrance  region  of  the  cell  is  conducted  to  the  exit 
location  where  it  is  removed  by  the  low  temperature  inlet  air. 


5.  List  of  symbols 


a  channel  thickness,  m 

R  ratio  gas-solid  heat-exchange  area/cell  area. 


<V 

</i, 

e 

i: 

F 

h 

/ 

k 

K 

n 

ne 

P 

/V  I 

<y 

R 

r 

K\ 

^pol 

S 

t 

T 

.v,  ’ 
V 


massic  specific  heat,  J  kg  1  K  " 1 
molar  specific  heat,  J  mol " 1  K  ~ 1 
hydraulic  diameter,  m 
emissivity,  - 
activation  energy,  J  mol  ' 

Faraday  constant,  A  s  mol '  1 
heat-transfer  coefficient,  VV  m  ’  K  1 
electrical  current  density,  A  m  ’ 
gas  heat  conductivity,  W  m " 1  K  1 
solid  heat  conductivity,  Wm'1  K  1 
molar  flow  rate,  mol  m  2  s  ~ 1 
number  of  electrons  transferred  in  the 
reaction,  - 
partial  pressure.  Pa 
reference  pressure.  Pa 
heat  flux,  W  m  2 
gas  constant,  J  mol  1  K1 
reaction  rate,  mol  m  *  s  1 
electrical  resistivity,  11  m 
polarization  resistance,  ft  nr 
cell  thickness,  m 
time,  k 

temperature.  K 
length  co-ordinates,  m 
potential,  V 


Greek  letters 


a 


$»«  Pi*  & 

A// 

AG 

AG° 


A 

v 

P 

a 


heat-transfer  coefficient  between  stack  and 
shell,  W  m“2  K~ 1 

coefficients  defined  in  Eq.  (9),  ft  K  m  , 
K 

thickness,  m 

enthalpy  change  of  the  over-all  reaction,  J 
mol"’ 

Gibbs  free  energy  change  of  the  over-all 
reaction,  J  mol  ” 1 

standard  Gibbs  free  energy  change  of  the 
over-all  reaction,  J  mol  ~  ‘ 
coefficients  defined  in  Eqs.  ( 10),  ( 1 1 ),  A 

m”2,  - 

stoichiometry  factor,  - 
stoichiometric  coefficient,  — 
density,  kg  m~3 

constant  of  Stefan-Boltzman,  W  m~ 2  K-4 


Subscripts 

an 

anode 

av 

average 

ca 

cathode 

g 

gas 

h 

layer 

/ 

chemical  species 

s 

solid 

shell 

shell 

.V 

.v-direction 

c-direction 

Superscripts 

oc 

open  circuit 

cc 

current  collectors 

ohm 

ohmic  effects 

pol 

polarization  effects 
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